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Protein denaturing induced by supercooling is interpreted as a process where 
some or all internal symmetries of the native protein are spontaneously broken. 
Hence, the free-energy potential corresponding to a folding-funnel landscape be¬ 
comes temperature-dependent and describes a phase transition. The idea that 
deformed vortices could be produced in the transition induced by temperature 
quenching, from native proteins to unfolded conformations is discussed in terms 
of the Zurek mechanism that implements the analogy between vortices, created 
in the laboratory at low energy, and the cosmic strings which are thought to have 
been left after symmetry breaking phase transitions in the early universe. An 
experiment is proposed to test the above idea which generalizes the cosmological 
analogy to also encompass biological systems and push a step ahead the view 
that protein folding is a biological equivalent of the big bang. 
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1 Introduction 


The folding of proteins can be regarded as a biological equivalent of the cosmo¬ 
logical big bang in that its end result is evident everywhere, in every living cell, 
but its beginnings and main cause are shrouded in mistery. This is of course 
mere analogy, but connections of the problem of protein folding with basic as¬ 
pects of fundamental physics are being discovered that call for a more serious 
interrelation. In particular, many proteins at sufficiently low temperature show 
a dynamical behaviour that matches that of glasses and spin glasses, while offer¬ 
ing a lot to be earned about the physics of complexity (Frauenfclder & Wolynes, 
1994). 

In this paper, we push forward the above analogy in still another important 
respect: we look at protein unfolding in supercooled pure solvent as being a 
process where some internal symmetries of the native protein are spontaneously 
broken. On the other hand, phase transitions induced by spontaneous symmetry 
breaking which are thought to have occurred in the early universe could have left 
behind long-lived topologically stable structures such as monopoles, strings, do¬ 
main walls and textures (Zcl’dovich et al, 1975; Kibble, 1976; Albrecht & Turok, 
1985; Vilenkin & Shellard, 1994). The possible role played by these structures in 
the generation of the present configuration of the observable universe has been 
very much discussed in recent decades (Zcl’dovich, 1980; Vilenkin, 1981; Turok, 
1983; Turok & Brandenberger, 1986). Therefore, the use of superfluid helium 
(Hendry et al, 1994; Buerle et al, 1996; Ruutu et al, 1996) or liquid crystals 
(Chuang et al, 1991; Bowick et al, 1994) to check in the laboratory the properties 
of the cosmological defects has become a matter of great interest. The possibil¬ 
ity of performing such cosmological experiments in helium was first suggested 
by Zurek (1985) who noted the analogy between cosmological strings and vortex 
lines in the superfluid. Here, this analogy is extended to also involve vortex lines 
filled with folded conformations of proteins rapidly denatured by lowering the 
temperature. 

2 Protein folding and spontaneous symmetry breaking 

Globular proteins which denature in water by lowering the temperature can be 
described as undergoing a phase transition driven by spontaneous breakdown of 
the internal symmetries hidden inside the core of their folded, native structure. 
Proteins in aquous medium interact with this through their hydrophobic and 
hydrophilic (polar) groups. At the sufficiently high temperatures for which no 
collective ’’coherent” effect can be expected to occur in the whole system, each 
protein molecule stabilizes by individually hiding most of its hydrophobic groups, 
placing them in the interior of the native, compact structure. The symmetries 
present in this structure can therefore be characterized by the number of contacts 
among its hydrophobic groups, and that number will tend to decrease as the 
protein unfolds (Dill & Su Chan, 1997). At lower temperatures, water can still act 
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like a collective field that attracts protein molecules to one another through the 
combined effect of the hydrophilic groups and water-to-water hydrogen bonding. 

The fundamental degrees of freedom here are the molecules themselves. A way 
to characterize the state with the thermodynamically most stable conformations 
is by the number density of protein molecules in the different conformations that 
exist, or that have ”condensed” into this state. Since at sufficiently high tem¬ 
peratures the unfolded protein molecules do not attractively interact with each 
other (i.e. their mutual weak attractive interaction is overcome by the combined 
effect of repulsive interactions among water and hydrophobic groups and thermal 
energy), the more unfolded conformations are formed the more energy it takes. 
However, because protein molecules are much less rigid in their unfolded confor¬ 
mations, water-mediated attractive interaction among the hydrophilic groups of 
two or more unfolded molecules is expected to be greater than that for protein 
molecules in the native conformation. Therefore, when the temperature becomes 
low enough to allow water molecules to tend to grouping around the polar groups 
of the protein, so that attractive interactions start dominating, the total energy of 
the system begins decreasing with the number density of unfolded conformations, 
until the system saturates in them. At such temperatures, one would expect all or 
some of the internal symmetries in the core of folded proteins to be spontaneously 
broken, creating a more stable ensemble of unfolded protein conformations. 

Let us consider the simplest form of a unidimensional section of the effective 
free-energy potential contribution satisfying the Anfinsen’s requirement (Anffii- 
sen, 1973) that native structures be thermodynamically stable states with con¬ 
formation at the global minimum of its accessible free energies when the inter¬ 
nal symmetry is not spontaneously broken, so as the new view that proteins 
change tertiary structure according to a folding funnel (Dill & Sun Chan, 1997) 
of their accessible energy landscape. In order to mathematically describe the 
phase change which is associated with protein unfolding induced by siipercooling 
the pure solvent, let us introduce the quantity T -the order parameter- which 
determines to what extend the geometrical distribution of hydrophobic groups 
in the unfolded phase at low temperature differs from that in the native phase, 
and adscribe a value T = 0 to this (symmetric) phase and a nonzero value of T 
to the (nonsymmetric) unfolded protein phase. The order parameter T and the 
free-energy potential defined in terms of T will respectively be associated with a 
given normal mode or combination of normal modes of the protein molecule, and 
with the corresponding section of the folding funnel landscape. 

Since the native protein is characterized by T = 0, near the low-temperature 
transition point one can write the free-energy potential as a power series in T. 
The number of relevant terms in this expansion should include those terms which 
are going to allow denaturation of the protein also at high temperature; i.e. we 
ought to cut off the expansion series only after the term T 8 to allow the existence 
of some local minima at 'f / 0 (visualizable at temperatures sufficiently higher 
than the transition temperature), while ensuring T = 0 to be the global minimum 
of the symmetric phase. Thus, a simplest free-energy potential that satisfies these 
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( 2 . 1 ) 


requirements can be written in the form 

1111 
= -«i^ 2 + -a 2 ^ 4 - -a 3 ^ 6 + -a 4 ^ 8 , 

2 4 o o 

where the aq’s are constants with the dimensions of a free energy, i = 1,...,4, 

such that a 3 > a 2 > a 3 > « 4 , and, since the states with 4/ = 0 and T ^ 0 

are assumed to possess different symmetry, the coefficient for the linear term 

has been neglected. Besides, we have restricted to the most interesting case 

which is prepared to accommodate transitions that follow a continuous process, 

rather than those occurring at isolated points, and therefore all the odd-power 

terms have been also neglected (otherwise, the transition points would be fixed 

by setting to zero more than one equation for the coefficients.) 

Apart of the global minimum at T = 0 , the potential free-energy ( 2 . 1 ) pos¬ 
sesses one local minimum and one local maximum. Spontaneous symmetry break¬ 
ing will occur in ( 2 . 1 ) when coefficient aq is allowed to take on negative values. 
The potential free-energy with coefficients aq constant, resulting from changing 
sign of aq in ( 2 . 1 ) will then show a local maximum at T = 0 , two local minima 
and one local maximum at T > 0. For the position of these extrema to be at 
real values of T in the temperature-independent approximation, the coefficients 
aq should moreover satisfy 


aq 


a% 


-aqaq ± I — a^aq 


a 2 ( a T 

+ 27 \“ 2 “ 4 “ T 


< 0 , 


( 2 . 2 ) 


where the lower signs correspond to potential ( 2 . 1 ) and the upper signs to the 
potential with negative aq. Manipulating these expressions, we obtain additional 
conditions on the aq’s: 


a 2 > 27aqa3, 


at > 


27 


-G2G4. 


(2. 3) 


If we consider the dependence of the free-energy potential on the thermody¬ 
namic quantities P, T of the protein, then the corresponding thermodynamic 
potential should be represented as a function V (P, T, T), where T must be deter¬ 
mined from the thermal equilibrium condition for which V is a global minimum. 
Restricting to processes taking place at constant pressure, it follows that aq(T) 
must vanish at the transition point because this term must be larger than zero 
in the symmetric phase (to have a global minimum at T = 0 ) and smaller than 
zero in the phase with broken symmetry to account for the existence of a global 
minimum at 4 / 0. Besides, for the transition point at which aq vanishes to 
be a stable state, it is necessary that the last term a^ be positive and constant. 
However, one can expect that the terms a 2 (F) and a 3 (T) still show a singularity 
where they change sign as function of T. These singularities would correspond to 
new transitions occurring at temperatures which must differ from the main tran¬ 
sition temperature. Thus, for a given P =const., one may write a 3 oc (T — T cl ), 
a 2 oc (T — T c2 ) and a 3 oc (T — T c3 ), and, given the shape of the potential, 
Pel < T c2 < T c 3 . 
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Once the temperature-independent potential (2.1) is fixed, the above results 
can also be obtained by using the following procedure. If we were dealing with 
a field theory in which coordinate was taken to be the matter field, then we 
had a precise prescription to investigate the symmetry behaviour of the effective 
theory with negative aq at nonzero temperature (Linde, 1979). However, what 
we, instead, have is a phenomenological theory with a potential linear in free 
energy which is not but a generalization of that of Ginzburg and Landau for 
superconductivity (Ginzburg & Landau, 1950). Field theory, in turn, should be 
expected to be nothing but a covariant generalization of the Ginzburg-Landau 
theory for a potential given in terms of energy densities rather than thermody¬ 
namic potentials (free energies). To obtain a temperature-dependent potential 
in our phenomenological theory using the field-theoretical prescription (Linde, 
1979), we note that the potential of the Ginzburg-Landau theory can be for¬ 
mally reproduced by regarding the order parameter as a field, and using then the 
field-theoretical prescription as applied to a ” fermionic-like” field, rather than the 
usual boson case. The ultimate reason to follow this formal procedure resides in 
the fact that the Ginzburg-Landau potential, or potential (2.1), is linear in free 
energy, rather than energy-density. 

Starting with the temperature-independent potential (2.1), we shall extend 
this trick to attain the corresponding temperature-dependent effective potential 
for our unidimensional symmetric slice of a protein folding funnel landscape. The 
Lagrange equation for the ’’field” coordinate T in the broken symmetry phase 
would be given by 

(□ + ctq - a 2 'b 2 + a 3 T 4 - a 4 T 6 ) T = 0. (2. 4 ) 

We then tentatively interpret the global minimum in (2.1) as defining a ’’con¬ 
formational vacuum”. Therefore, in Eqn. (2.4) we should first insert a shift 
T —> T = T + cr(T) to recover the picture where we can define usual creation 
and annihilation operators with vanishing conformational vacuum expectation 
for the ’’field” T, and then take the Gibbs average (Landau & Lifshitz, 1975) 

(...) = S ^j xp ^ , with H the Hamiltonian and T the temperature. We then 

obtain 

□<7(T) + a l( j{T) - a 2 (a{Tf + 3u(T)(T 2 )) 

+ a 3 (o-(T) 5 + 10u(T) 3 (T 2 )) - a 4 (u(T) 7 + 21a(T) 5 (T 2 >) = 0, (2. 5 ) 

where we have used the condition (T) = 0 for the new ’’conformational vacuum”, 
and restricted to work in the lowest order in a 2 , a 3 and a 4 , where the quantities 
(T 3 ), (T 5 ) and (T') vanish, and the contribution from (T 4 ) can be discarded. 

Since in the harmonic approximation at constant cr, ai('L 2 ) is linear in the 
free-energy, this quantity can now be calculated according to the fermionic rules 
as 

“■<*’> = (5 (2 - 6) 
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in which c u p is the energy of the particles with momentum p and mass m, and n p is 
the occupation number. Taking n p = (exp + 1 j , discarding the temparature- 
independent term — / d 3 p/ 2 cu 2 , which can be eliminated by mass renormalization 
at T = 0, and recalling (Linde, 1979) that all interesting effects should take place 
at T m, expression (2.6) becomes 


«i(d/ 2 ) 


1 r°° dp 

2n 2 Jo er + 1 


Tin 2 

’ 


(2. 7) 


so that we would finally have a Lagrange equation that again corresponds to an 
effective potential free-energy given by 


V = -ictieiT 2 + ^a 2 e 2 ^ 1 - ia^'L 6 + 
where the relative temperature ej are defined as: 


T ci - T 

J , 3 = 1,2,3 


e,- = m 
2 T ■ 

± CJ 


with the critical temperatures 


rp _ f rri _ r ^ 1^2 rj, rp _ r ^ 1^3 rp 

J cl ~ Jl-, ^c2 ~ J 2- > J clj J c3 — /3- > J c2, 

CC2 CC3 0:4 


040:2 


O1O3 


( 2 . 8 ) 

(2. 9) 


( 2 . 10 ) 


in which conditions (2.3) have been used and the dimensionless numerical coef¬ 
ficients /i, f- 2 and f 3 are all of order unity. The same approximate result would 
also be obtained if we took into account the discarded terms proportional to (T 4 ) 
in (2.5), since such terms, which can be calculated by the above procedure, would 
contribute T cl and T c2 with a factor of order unity because of conditions (2.3). 
These results can readily be generalized to the more realistic case of a potential 
with an arbitrary number of higher-order potential terms, the highest of which 
being even and positive, and terms with odd power of T which break invariance 
under T —» — T. For the purposes of this work, however, it will suffice working 
with the simplest potential slice ( 2 . 8 ). 

Note that (2.8) reduces in fact to the typical thermodynamical potential of 
the Ginzburg-Landau theory (Ginzburg & Landau, 1950) in the limit < 43 , aq —> 0, 
with <44 going to zero more rapidly than a 3 . The considered phase transition for 
protein folding can therefore be regarded as a generalization from the second- 
order phase transitions taking place in the phenomenological theory of supercon¬ 
ductivity. At T c i there will be a second-order phase transition from a potential 
with a shape as in ( 2 . 1 ) to the corresponding potential with the broken internal 
symmetries. According to (2.8), other transitions can also occur at the generally 
higher critical temperatures T c2 and T c3 , these taking place between states with 
the folded protein at the global minimum. 
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3 Cosmology in a protein 

Like in superconductors and superfluids (Tilley & Tilley, 1986), or field theories 
(Nielsen & Olesen, 1973), in the present phenomenological model strings filled 
with a protein in its symmetric, native phase will be able to form when the order 
parameter T is allowed to be complex 

T = |T|e* e , (3. 1 ) 

where 9 is a fixed phase. The possibility of formation of vortex lines in the con¬ 
sidered phase transition in proteins required therefore having a two-dimensional 
section of the inverted folding funnel (Dill & Su Chan, 1997) that corresponded 
to the breakdown of a complete given symmetry. Folded proteins contain a lot 
of internal symmetries which can be visualized as contacts between hydrophobic 
groups in the protein core (Dill, 1995; Dill & Su Chan, 1997). Each of these 
symmetries can be independently broken by variation of one, two or possibly 
more normal coordinates. In the first case, a domain wall could be formed; when 
variation of two normal coordinates were required to break down the hydropho¬ 
bic contacts defining one of such internal symmetries, then one would expect the 
formation of stringy topological defects, a case which we assume to hold hereafter. 

Generally, strings can be formed because in the broken-symmetry phase where 
T < T c i, the first term in (2.1) is negative and, since the last term of this potential 
is always positive, V has the shape of a wrinkled mexican sombrero; then, the 
manifold of the most stable conformations (conformational ” vacuum” manifold) 
is not simply connected and, therefore, possesses nontrivial loops. By the first 
homotopy group (Kibble, 1980; Vilenkin & Shellard, 1994) it follows that axisym- 
metric topological defects (vortex lines) should then be expected to form along 
the phase transition at temperature T c \. This can be most clearly understood by 
looking at the broken-symmetry phase with inverted folding funnel potential as 
a quantum system. 

Below the critical temperature T cl , a temperature-dependent fraction of pro¬ 
tein macromolecules become trapped in their various unfolded conformations and, 
together with the solvent molecules, will form a phase which can be described by a 
wave function T (r). The properties of the protein-solvent complex in the broken- 
symmetry phase can be understood assuming that this wave function satisfies a 
wave equation of the form 

(3 - 2) 

where m is taken to be the mass of the complex formed by a protein macro- 
molecule and its hydration shell, and yU is the chemical potential, i.e. the energy 
gained by the system when one molecule of protein is unfolded at constant volume 
and entropy. 

Let us assume that the free energy can be expanded in powers of |T| 2 (which 
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(3. 3) 


would now play the role of the order parameter) and has the form: 

F(r) = 0l (T)|>P| 2 + ia 2 (T)|'If - ia 3 (T)|1f + ia 4 U| 8 , 

where aq(T) = ejOLj- 

The chemical potential p can be evaluated by using (3.3) so that the Schrdinger 
equation becomes 

OiTf 12 

in-zr = ——v 2 T + aATW + a 2 (T) |T| 2 T - a 3 (T)|T| 4 T + a 4 |^| 6 ^. (3. 4 ) 
at 2 m 

Eq. (3.4) can now be re-scaled to assure a dimensionless form, which is given by 


if] = -V 2 r] + (-1 + |?/| 2 - |?/| 4 + |? 7 | 6 ) r], 


by introducing the quantities 


h 

«icni 


Tqi 

l e ll 
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h 

^2m\ai(T) \ 


Jo_ 

v^T’ 





2 CX\{T') OL\£\ 

V = ~a,MT) = (T)’ 


(3. 5) 

(3. 6) 
(3. 7) 
(3. 8) 


where T\ is to be interpreted as the relaxation time, £i as the correlation length, 
and the crj’s as equilibrium linear densities characterizing the broken-symmetry 
phase for the unfolded protein. In these definitions, is a time parameter that 
also characterizes the phase with broken symmetry. 

When we restrict ourselves to deal with the time-independent solutions to Eq. 
(3.5) in cylindrical coordinates (p,4>,z), i.e. 


dr] d 2 i] 
pdp dp 2 

one has the trivial solution 


\p | 2 = 1 and the axisymmetric solution 
V = f(p) exp (in0), 


(3. 9) 


(3. 10 ) 


with n a whole number. The radial part of (3.10) is regular for /) < 1, where 
f(p) ~ p n , and approaches the equilibrium value of the trivial solution at p 3> 
1, for which / 2 (p) ~ 1 — n 2 /2p 2 . Thus, solution (3.10) represents a stringy 
topological defect which, for n = 1, is still of the kind first considered by Ginzburg 
and Pitaevskii in the context of the phenomenological theory of superfluidity 
(Ginzburg & Pitaevskii, 1958). Around the axis of symmetry, there will be a 
flow caused by the phase gradient, with velocity v = TiSJO/m = h/mr. The 
resulting vortex will be filled with protein molecules keeping the native structure, 



or some partially folded conformations, and will be surrounded by a distribution 
of distinct unfolded configurations. It appear that a certain nonzero density of 
these vortex lines would be formed if the second-order phase transition is induced 
by a sufficiently rapid temperature quenching. 

However, the isotropy assumption for the order parameter, by which varia¬ 
tion of the accessible free energy respect to normal coordinates T i and T 2 are 
regarded to be equal, is obviously not but an approximation. Strictly speaking, 
the cylindric coordinates cannot be circular, but only closed on the ^-constant 
sections. This would lead to stringy solutions corresponding to a cylinder de¬ 
formed on its ^-constant sections and closed flow streamlines no longer perfectly 
circular. These deformed solutions would satisfy (3.10) only approximately. 

We shall estimate the density of vortex lines formed during the phase tran¬ 
sition of a protein induced on a dynamical time scale by a rapid decrease of 
temperature. In the vecinity of T c i, the relaxation time scale T\ (i.e. the time 
scale on which the order parameter can adjust to the new thermodynamic pa¬ 
rameters) will become quite longer than the time on which the quench proceeds. 
Then, the order parameter T will be practically frozen on that time scale and 
correlated over large domains where £1 —» oo (impulse regime). Sufficiently far 
from T cl , T\ will be much smaller that the quench time and T will become at 
an equilibrium configuration with determined by the instantaneous value of ei 
(adiabatic regime) (Zurek, 1996). 

If the transition is made fast enough, then the growing regions of the new 
phase met and coalesced, making it impossible to avoid the creation of topological 
defects. That will happen on the boundary between the impulse and adiabatic 
regimes, at the freeze out time t (Zurek, 1996), if 

Rb > L (3. 11 ) 


where is the effective radius of the region with the starting symmetric phase, 
and is the characteristic correlation length at the freeze out time t. 

At a sufficiently rapid pace, on a quench time scale tq (controlled by the rate 
at which the temperature is lowered), one can assume (Zurek, 1996) that is 
proportional to time, e\ = W, in the vecinity of T cl . Thus, at t there will be a 

freeze out temperature ej = £/tq, so that using (3.6), we obtain t = ^/Tol tq. The 
transition between the adiabatic and impulse regimes will therefore take place 
at the relative temperature ef = e(f) = \Jt 0 i/tq, corresponding to a freeze out 
correlation length which, using (3.7), can be expressed as 
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Jo_ 



\ 

Am 2 a i) 


1 

4 


(3. 12 ) 


Since an corresponds to an free energy which may be associated with a large 
number of the protein internal motions, such as unfolding of a-hclices (Creighton, 
1992) and hindered rotations, each with energies typically of the order 10 -13 erg., 
for most proteins will take on values within the interval (10~ 10 — ICC 13 ) erg., 
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and hence £1 will have values in the interval (10 6 — 10 7 )tq cm. For sufficiently 

rapid, but still not very fast temperature quench, we then obtain (i<l cm, and 
therefore there will be a copious initial production of vortex lines with density 
(Zurek, 1996). 

We propose finally an experiment to check the possible formation of vortex 
lines and other topological defects during the phase transition leading from folded 
to unfolded conformations in apomyoglobin (aMb). The transition will be induced 
by rapid supercooling, following excitation by focused pulses of a tunable infrared 
laser. The experiment is a variant of that recently carried out by Ballew et 
al. (1996) in order to follow the main fastest events of the aMb folding. It 
would consist of first denaturing aMb by supercooling it in water or D 2 0 at 
around -10 C in a perfectly insulated cell equipped with a high precision device 
to monitor variations of energy locally. This could be e.g. made by measuring the 
transmission of a calibrated low-power diode laser focused at the given region. 
The sample would be held in a cuvette with optical path length of the order 1 mm, 
cooled by thermoelectric devices keeping the temperature constant at around -10 
C, or inducing a fast recovery of that temperature, in the region previously heated 
by the infrared laser, with a thermistor feedback loop. Two counterpropagating 
beams from an infrared laser would here coherently and uniformly heat a cylindric 
region with section ~1 mm 2 , of the previously cooled sample, inducing the folding 
of aMb in the resulting irradiated volume. The pulsed infrared laser delivered 
short pulses (~ 1 ns) at wavelengths able to excite a near-infrared mode of water 
only (Ballew et al, 1996). After absorbing the infrared photons, water will almost 
instantly transform all absorbed energy into heat that locally induced naturation 
of aMb. 

After the laser excitation, the created interphase would inmcdiately start trav¬ 
eling inward the irradiated volume, quenching it to the denatured conformations 
in a time tq short enough to satisfy (3.11), so creating a tangle of topological de¬ 
fects filled with folded or partly folded aMb conformations. The presence of these 
defects would finally be detected by measuring the deficit in energy released as, 
in the irradiated volume, aMb molecules in the symmetric phase pass back into 
the low-temperature unfolded phase. By ascribing this deficit to the formation 
of vortices and other topological defects (Buerlc et al, 1996), one could infer the 
resulting defect density. 


Acknowledgements. P.F.G.-D. acknowledges DGICYT by support under Re¬ 
search Project No. PB94-0107. 
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